Comparison of Supernovae datasets Constraints on Dark Energy 
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Cosmological measurements suggest that our universe contains a dark energy component. In 
order to study the dark energy evolution, we constrain a parameterized dark energy equation of 
state w(z) = wo + WijXi us i n S the recent observational datasets: 157 Gold type la supernovae and 
the newly released 182 Gold type la supernovae by maximum likelihood method. It is found that 
the best fit w(z) crosses —1 in the past and the present best fit value of w(0) < —1 obtained from 
157 Gold type la supernovae. The crossing of —1 is not realized and wo — — 1 is not ruled out in 
lcr confidence level for the 182 Gold type la supernovae. We also find that the range of parameter 
wo is wide even in lcr confidence level and the best fit w(z) is sensitive to the prior of Q m . 

PACS numbers: 98.80.Es, 98.80.-k, 95.36.+X 



The universe is currently accelerating revealed by re- 
cent observations of type la supernovae (SNe Ia)[H, 0| 
and the acceleration has been attributed to a mys- 
terious component dubbed dark energy. Although a 
cosmological constant A is the simplest candidate for 
dark energy and appears to explain current observations 
satisfactorily [H, 0, it suffers from fine tuning and coinci- 
dence problems @- These difficulties have lead to a large 
variety of proposed models with time-depended dark en- 
ergy, such as quintessence B , phantomBl , quintom[l], K- 
essenceQ, tachyonic matter[l(| and so on. In all the 
above models, the evolution history of the universe and 
the nature of dark energy strongly depend on the models, 
such as the dynamics of the extra energy component and 
the gravity theory used to solve the cosmological prob- 
lem. However, there is another approach to study dark 
energy properties directly by observations in a model 
independent manner. In this method, the equation of 
state of dark energy w(z) = p(z)/p(z) could be tested 
by cosmological observations to understand the gravita- 
tional properties of dark energy and its evolution. The 
aim of this paper is to compare the 157 Gold SNe la 
(SN157)[H| with the latest 182 Gold SNe la (SN182)H2 
in constraining dark en erg y with a parameterized equa- 
tion of state O, H EE 



w(z) = wo + iv\ 



1 + 2 



(1) 



where w and w\ are constant, z is redshift. Since the 
spatially flat universe f2 m + VLde — 1 agrees well with 
observations [1,01, in this paper, we just consider the flat 
Friedmann- Robertson- Walker cosmology. 

In a flat universe with Eq.([T]), the Friedmann equation 
can be expressed as 

h\z) = H 2 [n m (i + zf + n DE (i + z )3(i+-o+-i) e ^Sf j 

(2) 

Then the knowledge of Q m and H(z) is sufficient to de- 
termine w(z). 

Supernovae observations provide the apparent magni- 
tude m(z) at peak brightness after implementing correc- 
tion for galactic extinction, K-correction and light curve 



width-luminosity correction. m(z) is related to the lumi- 
nosity distance cIl{z) through 

m th (z) = 51og 10 (I? L (z))+M + 51og 10 (^L-) + 25 

= 51og 10 (£> i (^)) + M + 42.38 -51og 10 /i (3) 

where Hq = hkm • s _1 • Mpc~ l , parameter M is the 
absolute magnitude which is assumed to be constant and 
Dl(z) is defined as 



D L (z) = H d L {z)/c 



= (1 + z) 



HqcIz 
o cH{z') 



(4) 



Supernovae data points are given in terms of the distance 
modulus 



fiobs(zi) = m obs (zi) - M 



(5) 



and theoretical model parameters are determined by min- 
imizing the quantity 



E (fJ.obs(zj) - HthiZj)) 2 
a 2 

i=l (obs-t) 



(6) 



where a 2 , .-, are the errors due to flux uncertainties, in- 

\obs;i) ' 

trinsic dispersion of SNe la absolute magnitude and pe- 
culiar velocity dispersion respectively. These errors are 
assumed to be gaussian and uncorrelated. The theoreti- 
cal distance modulus is defined as 

lHh{z%) = m th (zi) - M 

= 5 log 10 (£>i,(z)) +42.38- 5 log 10 ft (7) 

and fi b s is given by supernovae dataset. 

Our calculations to limit parameters are based on the 
Markov Chain Monte Carlo program CosmoMC[l7j. We 
modified the code to make equation of state w be redshift- 
depended instead of constant —1 and add parameter wq 
and w\ as slow parameters to the cosmological parame- 
ters space. Using the maximum likelihood method with 
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FIG. 1: The contours show 2-D marginalized la and 2a confi- 
dence limits in the (wi, w ), (wi, f2 m ), (Q,de, fim), (wo, f2m) 
plane from SN157. 
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FIG. 2: The contours show 2-D marginalized la and 2a confi- 
dence limits in the (vii, w a ), (tui, f2 m ), (Qde, f2,„), (w , f2 m ) 
plane from SN182. 




FIG. 3: The best fits of w(z) with la errors (shaded region). 
Q m is set to the best values obtained from SN157 or SN182. 




FIG. 4: The best fits of w(z) with la errors (shaded region) 
with a prior Q m — 0.27. 



TABLE 1. The minimum \ 2 an< A mean fit values from 
SN157 and SN182. 



data Xmin 


wo 




Wl 






SN157 89.947 


-1.856 ± 


0.519 


-0.015 ± 


0.269 


0.437 ± 0.073 


SN182 81.523 


-1.637 ± 


0.496 


-0.180 ± 


0.081 


0.435 ± 0.058 



SN157 and SN182, we obtained confidence level of equa- 
tion of state and matter density. The contour plots are 
show in FigQ] and Fig(2l the mean fit values, la con- 
straints on parameters wo, wi and fi m with standard 
deviations are list in Table 1. In Fig Q] we show results 
for the parameterized equation of state for the SN157 
dataset. We find the best fit value of (wi,wo,O m ) to 
be (0.335,-2.373,0.485) and see that the (wi,w ) panel 
has two island separately as well as the (wi,Q m ) panel 
which means that the SN157 dataset can not give a fine 
constraint on wi . Figf2] shows the results calculated by 
the SN182 dataset. The best fit value of (wi,wo,£l m ) is 
(-0.244,-1.728,0.462). 

Fig[3] shows la errors of the best fit w(z) calculated 
using the covariance matrix 18]. The best fit w(z) con- 
strained by SN157 can cross —1 in the evolution and in la 
confidence level, the LCDM parameter (wq = —l,Wi = 
0) do not coincide with the best fit value of w(z) at 
present. But the best fit dynamical w(z) obtained from 
SN182 does not cross —1 for the whole redshift range 
< z < 1.75, wo = — 1 is in la confidence contour and 
near the boundary. During our calculations, we also find 
that the best fit w(z) and its errors are sensitive to the 
value of Sl m . When we set r2 m = 0.27 as a prior instead 
of the best fit value mentioned before, both of the best fit 
w(z) can cross —1 and its la error region shrink sharply 
as showed in FigfJ] 

In summary, we have provided a comparative analysis 
on constraining parameters w\ , wq , Q m and the proper- 
ties of dark energy with two cosmological observations: 
SN157 and SN182. It shows that the SN157 dataset can 
not give a fine constraint on w\ and the best fit w(z) can 
cross —1 in the evolution, but the SN182 dataset gives 
contrary result. We also find that the range of parame- 
ter Wq is wide even in 1<t confidence level and the best fit 
w(z) is sensitive to the prior of Q m . In the future, we can 
combine other cosmological observations, such as Cosmic 
Microwave Background anisotropics [4j, the SDSS baryon 
acoustic peakyJl, the cluster baryon fraction (CBF),20] 
or the linear perturbations growth rate at z = 0.15 ob- 
tained from the 2dF galaxy redshift survey (PGR) 21] to 
constrain the matter density and the evolution of dark 
energy. We hope that the best fit values of cosmological 
parameters and the confidence level could be enhanced 
which will help us to understand the nature of dark en- 
ergy 
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